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ABSTRACT 

Based on the latest MLCS17 SNe la data provided by Hicken et al. (2009), 
together with the Baryon Acoustic OsciUation (BAO) and strong gravitational 
lenses (SGL), we investigate the dark energy equation-of-state parameter for both 
constant w and time-varying w{z) = Wq + Waz/{1 + z) in the flat universe, and 
its correlation with the matter density VLm and Hubble constant h. The con- 
straints from SNe data alone arrive at: (a) the best-fit results are {QM,w,h) = 
(0.358,-1.09,0.647), while both Qm and w are very sensitive to the difference 
Ah = h — hoi the Hubble constant deviating to the prior input h = 0.65; (b) the 
likelihoods of parameters are found to be: w = — O.SStgJg and ^Im = 0.3611q^I, 
which is consistent with the ACDM at 95% C.L.; (c) the two parameters in 
the time-varying case are found to be {wo,Wa) = (— O.TS^qJ^, 0.84^Jq^34) after 
marginalizing other parameters; (d) there is a clear degeneracy between constant 
w and Qm, which depresses the power of SNe la to constrain both of them; 
(e) the likelihood of parameter Wa has a high non-Gaussian distribution; (f) 
an extra restriction on Qm is necessary to improve the constraint of the SNe 
la data on {wq, Wa)- A joint analysis of SNe la data and BAO is made to 
break the degeneracy between w and Qm, and it provides a stringent constrain 
with the likelihoods: w = -0.88+[J:[J^ and Qm = 0.29;[}:|53. For the time- varying 
w{z), it leads to the interesting maximum likelihoods Wq = —0.94 and Wa = 0. 
When marginalizing the parameters Qm and h, the fitting results are found to 
be {wo,Wa) = (— 0.95^o;^g, 0.41]'^o;gg). After adding the splitting angle statistic 
of SGL data, a consistent constraint is obtained (fi^, w) = (0.298, —0.907) and 
{wo,Wa) is further improved to be {wo,Wa) = (— 0.92]'^q;Jq, 0.35]'^q;54), which indi- 
cates that the phantom type models are disfavored. 
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1. INTRODUCTION 

Analysis of the distance modulus versus redshift relation of type la supernova (SNe la) 
provides a direct evidence that the universe expansion is accelerating in the last few billion 
years(e.g., Riess et al. 1998, 2004, 2007; Wood-Vasey et al. 2007; Kowalski et al. 2008; 
Hicken et al. 2009; Biswas and Wandeltl 2009). This cosmic image is also supported by many 
other cosmological observations, hke the Cosmic Microwave Background(CMB)(Hinshaw et 
al. 2009; Komatsu et al. 2009), the Baryon Acoustic Oscillation (BAO) measurement 
(York et al. 2000; Eisenstein et al. 2005) and the weak gravitational lenses (Weinberg 
and Kamionkowski 2002; Zhan and Knox 2006). Based on the Friedmann equation, the 
acceleration can be explained through introducing a negative pressure component in the 
universe, named dark energy, which is nearly spatially uniform distribution and contributes 
about 2/3 critical density of universe today. To reveal the property of dark energy, the 
most of studies, either theoretical models or experiment data analysis, are focused on its 
equation-of-state parameter w — p/p(Albrecht et al. 2006). Here we shall utilize the latest 
SNe la data provided by Hicken et al. (2009) with using MLCS17 hght curve fitter, together 
with the sphtting angle statistic of strong gravitational lenses(SGL; Zhang et al. 2009) and 
the baryonic acoustic oscillations (BAO; Eisenstein et al. 2005) to investigate the constraint 
for the parameter w of dark energy in the flat cosmology. The influences of Hubble constant 
h and matter density Qm on the fitting results are carefully demonstrated. 

By far, all observed data are consistent with the ACDM cosmology, with dark energy 
in the form of a cosmological constant A. However, this model raises theoretical problems 
related to the fine tuned value (see e.g. Padmanabhan 2003). Many other theoretical models, 
like quintessence models and phantom model(Ratra and Peebles 1988; Caldwell, Dave, and 
Steinhardt 1998; Caldwell 2002) , reveal that dark energy might be a dynamical component 
and evolves with time. It is usual to parametrize dark energy as an ideal liquid with its 
equation-of-state(EOS) parameter w{z) = wo + Wa z/{l + z) (ChevaUier et al. 2001; Linder 
2003). Conveniently, it includes the case of a constant EOS with (wq — w,Wa — 0), and 
the ACDM model {wq — —l,Wa — 0). Then theoretical models can be classified in a phase 
diagram on the {wo,Wa) plane (see e.g. Barger et al. 2006; Biswas and Wandelt 2009). 
Thus the accurate measurement of the parameters (wo-Wa) is helpful for testing a certain 
theoretical models. The current allowed regions of {wo,Wa) given by the most observations 
or their combinations remain surrounding the crucial point {wq = —l,Wa = 0), which is 
the common point in the phase diagram of different classified models. Therefore, the final 
judgment of models can not be made and more careful works are still needed. 

In practice, the measurements of other cosmological parameters, like Hubble constant 
h and matter density parameter Qm, shall affect the determination of the EOS parameter w 
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of dark energy. Hubble constant h gives the rate of recession of distant galaxies and scales 
the expansion of the present universe. Parameter Qm describes the present matter density 
relative to the critical density, which is essential to understand the structure formation. In the 
flat cosmology, Qm is directly related to the present density of dark energy hj Q^e = 1 — 
The precise measurements on the two cosmological parameters are still processing. The 
results from WMAP five-year observations are h = 0.719 ± 0.03 and Q = 0.258 ± 0.03 
in the ACDM cosmology, respectively (Hinshaw et al. 2009). The SNe la also have been 
used to determine the Hubble constant, calibrated with their peak luminosities or Cepheid 
variables in nearby galaxies (Schaefer 1996; Sandage, et al. 1996; Branch 1998; Gibson 2001; 
Sandage et al. 2006; Riess et al. 2009). Because the Hubble constant degenerates with 
the distance modulus of SNc la. the correct value of the absolute magnitude of the fiducial 
SN la or Hubble constant is not relevant for the dark energy analysis which only make 
use of differences between SNc la magnitudes(Riess et al. 2004). In this case, the value 
h = 0.65 is often used as a prior input(Hicken et al. 2009). However, we can still investigate 
the sensitivity of the analysis results on the change of the Hubble constant, which is not 
presented in the literatures. 

The SNe la has homogeneity and extremely high intrinsic luminosity of peak magnitude 
and thus is widely used to measure the cosmological parameters(e.g. Riess et al. 2004; 
Wood-Vasey et al. 2007; Kowalski et al. 2008; Hicken et al. 2009). With given density of 
dark energy in the universe today p{z = 0), the change of equation-of-state parameter w will 
bring the change of its density p{z) at the redshift z and then the distance d{z). Inversely, 
measurement of the redshift z of SNe la and corresponding distance d{z) can constrain 
the dark energy. In spite of the high accuracy of SNe la measurement, its potential of 
constraining dark energy is not very strong due to the degeneracy of w and CIm- Therefore, 
a combining analysis with other observations is often useful. 

We will also use the summary parameters of the baryonic acoustic oscillations as reported 
in previous studies(Eisenstein et al. 2005). The large-scale correlation function of a large 
sample of luminous red galaxies has been measured in the Sloan Digital Sky Survey and a 
well-detected acoustic peak was found to provide a standard ruler by which the absolute 
distance of z = 0.35 can be determined with 5% accuracy, independent on the Hubble 
constant h. This ruler is a prior and can be used to constrain dark energy(e.g. Porciani 
and Madau 2000; Huterer and Ma 2004; Chae, 2007). We shall also use strong gravitational 
lensing statistic observation, which provides us a useful probe of dark energy of the universe. 
Mainly through the comoving number density of dark halos described by Press- Schechter 
theory and the background cosmological line element, dark energy affects the efficiency with 
which dark-matter concentrations produce strong lensing signals. Then through comparing 
the observed number of lenses with the theoretical expected result as a function of image 
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separation and cosmological parameters, it enables us to determine the allowed range of the 
parameter w. The constraint process also depends on the density profile of dark halos. Here 
we shall use the two model combined mechanism to reproduce the observed curve of lensing 
probability to the image splitting angle (Sarbu, Rusin and Ma 2001; Li and Ostriker 2002; 
Zhang et al. 2009). The redshift of CMB is above 1000 and far larger than 1, and there is no 
other observation to fill up this redshift gap, thus we would not adopt the CMB data in the 
present analysis and limit our study on dark energy to the redshift region of 2; ~ 1, which is 
the characteristic redshift scale of SNe la, BAO and SGL statistic. 

In our recent work(Zhang et al. 2009), we have present the constraint on the dark 
energy from the SGL splitting angle statistic. In this paper, by taking the latest analyzed 
SNe la data(Hicken et al. 2009), the baryonic acoustic oscillations (Eisenstein et al. 2005) 
and the CLASS statistical sample(Browne et al.2003), we shall make an joint analysis to 
constrain the dark energy equation of state parameter w and the matter density JIm, and to 
study the possible influence of the Hubble constant h. Two model independent assumptions 
of dark energy are considered: constant w and the time- varying parameterization w{z) = 
Wo+Wa zj (1+z). We mainly highlight two issues which have not previously been illuminated. 
First, we carefully study, based on the latest SNe la data, the constraints for the dark energy 
EOS w{z) and the influences of the matter density and Hubble constant h. Second, we 
investigate the joint analysis of SNe la data, BAO and SGL statistic in detail. Our paper 
is organized as follows: Sect. 2 shows the constraint by the latest SNe la data on dark 
energy. Sect. 3 describes the joint analysis of the SNe la, the BAO and the SGL statistic, 
more stringent constraints on dark energy and matter density are resulted, and the influence 
of Hubble constant is explicitly demonstrated. The conclusions are presented in the last 
section. 



2. DARK ENERGY CONSTRAINTS BY THE LATEST SNe la DATA 

As the standard candles of the cosmology, the SNe la is used to study the geometry and 
dynamics of the universe with redshift z < 1.7. In determinations of cosmological parameters 
about the accelerating expand and dark energy, the SNe la remains a key ingredient. In 1998, 
the SNe la measurement provided the flrst direct evidence for the presence of dark energy 
with the negative pressure. Then many SN la observations have been done and the total 
number of SNe la sample increases quickly. The SN la compilations arc often consist of 
high-redshift {z ~ 0.5) data set and low-redshift {z ~ 0.05) sample at the same time (e.g. 
Riess et al. 1998; Perlmutter et al. 1999; Wood-Vasey et al. 2007). When combining several 
independent group's SNe la data sets into one compilation, the consistent analysis method 
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of light curves and the selection of supernova are crucial. For a certain sample, the different 
light curve fitter and corresponding different selection of supernova can lead to different 
constraints on the cosmological parameters(e.g. Hicken et al. 2009). 

The fitting results of cosmological parameters from different SN la compilations show 
an obvious difference which could be huge in some cases. For example, the sample provided 

by Wood-Vascy ct al. (2007) led a cosmology with near zero matter component VLm ~ 0, 
while the best-fit result of VLm by Ricss ct al. (2004) was found to be ~ 0.5. The fitting 
constant w of dark energy for the former sample was found to be to ~ —0.75 and for the latter 
sample, it was found to be w ~ —2.0. On the other hand, it is obvious that larger sample of 
SNe la is more powerful for constraining cosmological parameters. Therefore it is useful and 
attractive to adopt the SNe la data as many as possible. There is a conventional method 
to combine several group's SNe la compilations, namely, by introducing an extra nuisance 
parameter in the statistic of every used SNe la sample and marginalizing them over in 
the fit, all statistics of samples can be summed into one total statistics (see e.g. Barger 
et al. 2006). The nuisance parameters are considered as analysis-dependent global unknown 
constants in the distances and are degenerate with the Hubble constant h. Although this 
combined mechanism is wildly adopted, the so-called analysis-dependent unknown constant 
is just an averaged effect of analysis-dependent uncertainties. 

For a certain SNc la, its 'measuring results' from independent groups often have visible 
differences from disparate light curve fitting functions and analysis procedures. Ricss ct al. 
(2004) and Wood- Vascy ct al. (2007) provided 157 golden sample and 162 SNc la sample, 
respectively. The two samples both adopt MLCS2k2 light-curve fitter and have 39 common 
SNe la, every SNe la of which has two sets of different redshift z and distance moduh ji. 
From the 39 common data, we can compare the two data sample directly and then obtain 
the distance uncertainty Oo^fj, — 0.21 and redshift uncertainty Go^z — 0.00029. Therefore it is 
clear that the analysis-dependent uncertainties from different group are quite large. Partially 
to solve this problem, Kowalski et al. (2008) provided the Union data set, a compilation of 
307 SNe la discovered in different surveys. The heterogeneous nature of the data set have 
been reflected and all SNe la sample are analyzed with the same analysis procedure. 

In the Union data set, all SNe la light curve are fitted by using the spectral-template- 
based fit method of Guy et al. (2005) (also known as SALT). There are other light curve 
fitters used in literatures, such as SALT2 (Guy et al. 2007), MLCS2k2 (Jha, Riess, and 
Kirshner 2007) with Ry = 3.1 (MLCS31) and MLCS2k2 with Ry = 1.7 (MLCS17). Hicken 
et al. (2009) compared these light curve fitters and found that SALT produces high-redshift 
Hubble residuals with systematic trends versus color and larger scatter than MLCS2k2, and 
MLCS31 overestimates host-galaxy extinction while MLCS17 does not. For a certain SNe 
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la, the analysis outcomes of different light curve fitters are not equal. Here we choose the 
SNe la compilation provided by using MLCS17 light curve fitter with the best cuts Ay < 0.5 
and A < 0.7 to constrain the dark energy. 

In the flat universe, the Friedmann equation are given by 

J (l-fiM)(l + ^)^^^+"'^ for constant w, 

^^^^^ ~ j (l-r]M)(l + 2)^^'+""+"'"^e-3"'«^/(i+^) for w{z)=wo + Wajlr,, 

with Hubble constant Hq = 100 h km s~^ Mpc^^. The influence of cosmological parameter 
w is focused on the dark energy density Qde{z) and then the Luminosity distance d^, which 
is defined as 

Analysis of the distance modulus versus redshift relation of SNe la can give us the 
information about the cosmological parameters. Distance estimates of SNe la are derived 
from the luminosity distance, di = {C/ A-nJ^Y^'^ where C and JF are the intrinsic luminosity 
and observed flux of the SNe la, respectively. It is usual to introduce the apparent magnitude 
m and absolute magnitude M. The apparent magnitude m of a celestial body is a measure 
of its brightness as seen by an observer on Earth and it is defined as m = log 5/jqq(^/JF°) = 
— 2.5 logiQ(jF/jF''), where is a reference flux, which is the zero point by definition (used to 
be the observed flux of Vega star). The absolute magnitude M of a celestial body outside of 
the solar system is defined to be the apparent magnitude it would have if it were 10 parsecs 
away, namely 



M = m(lOpc) 



2 



-2.51ogio[^(10pc)/^°] = -2.51ogio[^ (M^^°'°] 



m — 2.5 log 



4 



25 (3) 



_(Mpc)2_ 

From the definition of the distance moduli fi = m — M, we have 

= 51ogrfL/Mpc + 25. (4) 

Using Equations ([T]), ([2]) and (jlj), we can relate the parameter w with the measured redshift 
z and distance moduli ^{z) of SNe la data. Define c/l = HodL, then we get 

/i + 5 log //o = 5 log dL + 25. (5) 
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The right side of the above formula doesn't depend on the Hubble constant and h directly 
degenerates with the distance modula /i. For the dark energy analysis, the SNe la data often 
use the certain value of h = 0.65 as the prior input and only make use of the differences 
of the logarithmic SNe la distances. However, we can still investigate the influence on the 
statistic analytic results when changing the Hubble constant h. Here Ah = h — h = /i— 0.65 is 
equivalent to the difference between the input Hubble constant and the real Hubble constant 
in the universe. Therefore although the certain value of the parameter h loses its physics 
meaning. Ah is still a physics quantity. The likelihood for the parameters {Qm,w, h) can be 
determined from a statistic 



where subscript i denotes the ith SNe la data and ai is the observed uncertainty. and 
/i"^ are the observed and theoretical distance moduli, respectively. 

Let us first discuss the constraints for the constant w case. Using the Powell mini- 
mization method(Press, et al 1992), we minimize the likelihood function of the three pa- 
rameters {QM,'w,h) and find that the coordinate of the best-fit point is {QM,w,h) = 
(0.358,-1.09,0.647). The fitted result of the Hubble constant h is very close to the prior 
input value h = 0.65 of Hicken et al. (2009) . It is noted that the value of Q^j is large, in 
comparison with Qm = 0.258 of the concordance cosmology provided by WMAP five year 
data(Komatsu et al. 2009). 

To see the correlations among the three parameters {QM,w,h) and their influences 
on the fitting results, we consider three different combination fittings which are shown in 
figure dJ From top to bottom in figure [H we show the fitting results of {flM,w), {w,h) 
and {^M, h) with other parameter being varying in the resonable range. In every panel, 
three pairs of contours are presented with a given cosmological parameter: for the upper 
panel, the deviation of the Hubble constant to the prior input Ah = h — 0.65 is taken to be: 
Ah = -0.03, Ah = -0.005 and Ah = 0.02, which are corresponding to h = 0.620, 0.645, 
and 0.670, from top to bottom respectively; for the middle panel, the matter density Qm is 
taken to be: Qm = 0.35, 0.30, and 0.25 from left to right, respectively; for the lower panel, 
the constant w of dark energy is taken to be: w = —0.7, —0.9, and —1.1 from left to right, 
respectively. The crosshairs in every panel mark the best fit points: for the upper panel, 
{nM,w) = (0.01, -0.41), (0.28, -0.87) and (0.48, -2.34) for h = 0.620, 0.645, and 0.670; for 
the middle panel, {w,h) = (-1.05,0.647), (-0.91,0.646) and (-0.81,0.645) for = 0.35, 
0.30, and 0.25; for the lower panel, {QM,h) = (0.19,0.644), (0.29,0.646) and (0.36,0.647) 
for w = —0.7, —0.9, and —1.1. It can be seen that the fitting results for (QmjW) from the 
SNe la data is very sensitive to the value of Hubble constant h: when h increases, the fitting 




(6) 



- 8 - 



Qm^w) increases(decreases) very rapidly. 

To investigate the constraint on the interested parameters, it is usual to marginalize 
other parameters. Figure [2] shows the {Qm,ui), {w,h) and {QM,h) contours from top to 
bottom, respectively. The unshown parameter in every panel has been marginalized through 
integrating it by applying the likelihood function L = exp(— The best-fit points are 
{nM,w) = (0.30,-0.92), {w,h) = (-0.86,0.645) and (^^M, /?-) = (0.35,0.647). The (1^m,w) 
contours have been shown in figure 4 of paper by Hicken et al. 2009, and our corresponding 
results are almost identical with theirs. Here we shall carefully discuss the information from 
the figures: (a) there is an obvious degeneracy of constant w and Qm, which depresses the 
power of SNe la to constrain both of them; (b) the fitting results are consistent with the 
ACDM cosmology at 95% C.L.. Comparing with the results obtained from the Union data 
set (Kowalski et al. 2008), we notice that the confidence regions in our present consideration 
are similar to those of the Union data set. Figure [3] shows the likelihoods of parameter Qm, 
w and h, in which the maximum likelihood points are located at Qm = 0.36, w = —0.88 and 
h = 0.646, respectively. It is interesting to notice that the parameter w is restricted to be 
from —2.0 to —0.5 and Qm is less than 0.5. Again, our fitted Hubble constant h is consistent 
with the prior input value h = 0.65 and its allowed region is very narrow. 

We now focus on the time-varying model w{z) = wq + Waz/{1 + z). Figure H] shows 
the contours of (wo,Wa), the best-fit point is {wq^Wo) = (—0.73,0.84). It is seen that the 
SNe la data alone have a poor constraint power on the parameter Wa- In figure we 
plot the contours of two parameters {flM,Wo), the best-fit point is found to be {flM,u!o) = 
(0.45, —0.68). It is shown that when Qm increases from 0.3 to 0.45, the allowed region of 
parameter wq is enlarged quickly. For a smaller Qm < 0.3, it leads to a much better constraint 
for the parameter Wq: Wq ~ —1.4 ~ —0.6. Figure [6] gives the contours of {QM,Wa), the best- 
fit point is {QMjWa) = (0.44, —4.63). It is noticed that when Qm increases from 0.34 to 0.5, 
the allowed region for the parameter Wa is enlarged rapidly. For a smaller Qm < 0.34, it 
also leads to a much better constraint for the parameter Wa- Wa ~ —3.0 ~ 2.5. From the 
figure [5] and figure [6l it indicates that an extra restriction on Qm is necessary to improve the 
constraint of the SNe la data on the parameters Wq and Wa- Figure [7] shows the likelihoods 
of parameters Wq and Wa, in which the maximum likelihood points are located at Wq = —0.8 
and Wa = 0.4, respectively. It can be seen that the parameter wq is limited in the region 
—2.5 < Wo < 0.5 and the likelihood of parameter Wa has a high non-Gaussian distribution. 
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3. JOINT ANALYSIS OF SNe la DATA, BAO AND SGL STATISTIC 

As we have shown in the previous chapter, over 300 SNe la observed so far are not 
sufficient for determining the cosmological parameters, especially for Wo and Wa- Many 
surveys(e.g. the Dark Energy Survey and Pan-STARRS) are proposed to obtain the SNe la 
sample with enlarged number and improved precision. Here we are going to constrain the 
dark energy through the combination of the SNe la data, the baryon acoustic oscillations as 
well as the SGL splitting angel statistic. 



3.1. Bciryon Acoustic Oscillations 

In the relativistic plasma of the early universe, ionized hydrogens (protons and electrons) 
are coupled with energetic photons by Thomson scattering. The plasma density is uniform 
except for the primordial cosmological perturbations. Driven by high pressure, the plasma 
fluctuations spread outward at over half the speed of hght. After about 10^ years, the 
universe has cooled enough and the protons capture the electrons to form neutral Hydrogen. 
This decouples the photons from the baryons, which dramatically decreases the sound speed 
and effectively ends the sound wave propagation. Because the universe has a significant 
fraction of baryons, these baryon acoustic oscillations leave their imprint on very large scale 
structures (about lOOMpc) of the Universe. 

The measurement of baryon acoustic oscillations was ffist processed by the Sloan Digital 
Sky Survey (SDSS; York et al. 2000) and Eisenstein et al. (2005) studied the large-scale 
correlation function of its sample, which is composed of 46, 748 luminous red galaxies over 
3816 square degrees and in the redshift range 0.16 to 0.47. The typical redshift of the sample 
is at 2; = 0.35. The large-scale correlation function is a combination of the correlations 
measured in the radial (redshift space) and the transverse (angular space) direction (Davis 
et al. 2007). Thus, the relevant distance measure is modeled by the so-called dilation scale, 
Dv{z) = [D\(z)z/H{z)Y^^, with comoving angular diameter distance Da{z) — dz' / H{z'). 
The dimensionless combination A(z) = Dv{z)^^/^ImH^/ z has no dependence on the Hubble 
constant h and is found to be well constrained by the SDSS data at ^ = 0.35. A standard 
ruler is provided as (Eisenstein et al. 2005) 



[H{zr)/Ho]y^ 



1 r 

Zl Jo 



dz 



H{z)/Ho 



2/3 



= 0.469 ±0.017 , (7) 



where zi = 0.35. This ruler is a prior and can be used to constrain dark energy. Then 
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the statistic is given by 



X i^M: W) = (8) 



3.2. SGL Splitting Angle Statistic 



The CLASS statistical sample has provided a well-defined statistical sample with N — 
8958 sources. Totally Ni — 13 multiple image gravitational lenses have been discovered and 
all have image separations < 3" (Browne et al. 2003). The SGL statistics arc sensitive 
to the equation-of-state parameter w of dark energy, which influences the number density 
of lens galaxies and the distances between the sources and lens. The probability with image 
separations larger than for a source at rcdshift Zg on account of the galaxies distribution 
from the source to the observer can be obtained by (Schneider et al. 1992) 



P(> A^) = 



dz 



[1 + zfn{M, z)a{> A9) dMdz , 



(9) 



where M is the mass of a dark halo, -Dl is the proper distance from the observer to the lens, 
n(M, z) is the comoving number density of dark halos virialized by redshift z with mass M 
and a is the cross section for two images with a splitting angle > A^. 

According to the Press-Schechter theory, the comoving number density with mass in the 
range (M, M + dM) is given by 



n{M,z)dM^^f{M,z)dM. 



(10) 



with the matter density of universe today po = ^MpcrUfi and the critical matter density at 
present Pcrit,o = SHq/^SttG). f{M,z) is the Press-Schechter function, and we shall utilize 
the modifled form by Sheth and Tormcn (1999) 



f{M,z)^ 



0.383 5r dA 



1 + 



A^ 



Vtt A'^dM 

dk 

A\M,z) = —A,{k,z)W'{k 



0.707(52 



0.3 



X exp 



0.707 fS, 



(11) 
(12) 



where A is the variance of the mass fluctuations (Eisenstein and Hu 1999) and parameter 
Sc{z) is the linear overdensity threshold for a spherical collapse(Wang and Steinhardt 1998; 
Weinberg and Kamionkowski 2002). 

For different density proflle of dark halo, the lensing cross section a can be calculated 
out based on the lensing equation. We shall use the combined mechanism of SIS and NFW 
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model to explain the whole experimental curve of strong gravitational lensing statistic. For 
that a new model parameter was introduced by Li and Ostriker (2002): lenses with mass 
M < Mc have the SIS profile, while lenses with mass M > have the NFW profile. Then 
the differential probability is given by 

dP/dM = dPsis/dMd{M^ - M) + dPNFw/dM§{M - M^) 

where d is the step function, ^{x — y) = 1, if x > y and otherwise. As the splitting angle 
A6 is directly proportional to the mass M of lens halos, the contribution to large A6 of SIS 
profile is depressed by M^. The lens data require a mass threshold ~ 10^^h~^MQ, which 
is consistent with the halo mass whose cooling time equals the age of the universe today. 

The likelihood function of the SGL splitting angle statistic is defined as 

Uw) = {l-p{w)f-''^llq.{w). (13) 

i=l 

p{w) and qiiw) represent the model- predicted lensing probabilities and the differential lensing 
probabilities, respectively. They are related to P in Equation ^ by an integration 

p{w) = Pobs(> ^^) = j j B "^^^^^^^^ i^{zs)dzdz, , (14) 

and 

cliw) = =J JB -^^^ ^{zs)dzdz, . (15) 

B is the magnification bias and can be found in our previous work(Zhang et al. 2009). i~p{za) 
is the redshift distribution of sources. Here we take the Gaussian model by directly fitting 
the redshift distribution of the subsample of CLASS statistical sample provided by Marlow 
et al. (2000), which is given by (Zhang et al. 2009) 



9{.Zs) = ^E=rexp 
v ZnA 



with A^, = 1.6125; a = 0.4224: A = 1.3761. 



3.3. Joint Analysis and Numerical Results 



In this section, we will investigate the constraint on the cosmological parameters from 
the joint analysis of (SNe + BAO) and (SNe + BAO + SGL), respectively. For the two(or 
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three) independent observations, the hkehhood function of a joint analysis is just given by 

L = i^SNe X -^BAO (X-^SGl) 

= exp(-xiNe/2) X exp(-x|Ao/2) (xLsgl). (17) 

The statistic significance Xbao xInb can be obtained by using Equations ([8]) and IQ, 
respectively. Lsgl is the likelihood function of SGL statistic and can be obtained by using 
Equation fll3p . The parameter h in the SNe la data and SGL data has different meaning 
and For the joint analysis of (SNe + SGL), the parameter h should be integrated separately: 
for the former, we integrate h in the range around 0.65 and for the latter, the integral range 
is 0.4 to 0.9. The BAO analysis has no dependence on the Hubble constant h. 

Let us first discuss the constraints for the constant w case from the joint analysis of 
(SNe + BAO). After the Powell minimization(Press, et al 1992), we get the best fit results 
of the three parameters {flM,w,h) = (0.287,-0.885,0.646). The fitted Qm decreases in 
comparison with the result of SNe data alone, and at the same time the fitted w increases 
due to their degeneracy relation. Then in figure [8], we show the constraints on Qm and the 
constant w from the joint analysis of (SNe + BAO). The parameter h has been marginalized. 
The best-fit results are {flM,'w) = (0.28, —0.88). The 95% C.L. allowed regions of constant 
w and Qm are found to be: —1.03 < w < —0.75 and 0.24 < Qm < 0.33, respectively. It 
is seen that: (a) the most allowed region of parameter w is above the line w = —1; (b) the 
fitting results are consistent with the ACDM cosmology at 95% C.L.. Comparing with the 
results of SNe la data alone, we see that the allowed regions for Qm and w are much reduced 
and the degeneracy between them disappears. 

In figure [H we show the constraints on Qm and the constant w from the joint analysis 
of (SNe + BAO + SGL). The parameter h has been marginalized again. The best fit result 
is {flM,w) = (0.29, —0.91). The 95% C.L. allowed regions of constant w and Qm are found 
to be: -1.06 <w< -0.77 and 0.25 < < 0.34. Comparing with the results of (SNe + 
BAO) case, it is seen that the results have only slight differences and the fitted w is found 
to be slightly smaller after adding the SGL data. 

Figure [TU] plots the likelihoods for the parameters Qm, w and h from the joint analysis 
of (SNe + BAO) and for the parameters Qm and w from (SNe + BAO + SGL), respectively. 
The maximum likelihood points are located at Qm = 0.29, w = —0.88 and h = 0.65 for 
(SNe + BAO) and = 0.296 and w = -0.91 for (SNe + BAO + SGL). It is interesting 
to find that the parameters w and Qm are restricted to the range: —1.17 < w < —0.67 and 
0.23 < Qm < 0.37. 

After marginalizing the cosmological parameters {QM,h), we obtain the constraint on 
{wo,Wa) in figure [m from the joint analysis of (SNe + BAO) and (SNe + BAO + SGL), 
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respectively. The crosshairs mark the best-fit point (wq, Wa) = (—0.95, 0.41) for the (SNe + 
BAO) case and (wo,wJ = (-0.92,0.35) for the (SNe + BAO + SGL) case. For the (SNe 
+ BAO) case, the 95% C.L. allowed regions for the parameters Wq and Wa are found to be: 
-1.22 < u;o < -0.66 and -0.92 < Wa < 1.59. For the (SNe + BAO + SGL) case, the 95% 
C.L. allowed regions for the parameters wq and Wa are found to be: —1.10 < wq < —0.72 
and —0.55 < Wa < 1.32. After adding the SGL data, the constraint on the parameter wq is 
improved moderately, but for the parameter Wa, the allowed region decreases by near half. 
The extra constraint power on the time- varying w{z) obtained through adding SGL data is 
due to the larger redshift < 2; < 3.0 of the galaxies in CLASS observational sample, in 
comparison with the reshift range of SNe data < z < 1.5 and the redshift of BAO z = 0.35. 
It can be seen for the both cases that: (a) the most allowed region of Wa is above Wa = 0; 
(b) in comparison with the cosmological constant {wo,Wa) = (—1.0,0.0), the joint analysis 
for both cases favors more positive {wo,Wa); (c) in comparison with the results of SNe la 
data alone, the constraint on Wa is much improved and wo also gets better constrained. 

Figure [12] plots the likelihoods of parameters wq and Wa from the joint analysis of (SNe 
-I- BAO) and (SNe + BAO + SGL), respectively. For (SNe + BAO) case, the maximum 
likelihood points are located at Wq = —0.94 and Wa = 0. Note that Wa = implies a constant 
equation-of-state of dark energy. For (SNe + BAO + SGL) case, the maximum likelihood 
points are found to be wq = —0.91 and Wa = 0.34. We see that the parameters wq and Wa 
are restricted to be: —1.20 < wq < —0.67 and —1.0 < Wa < 2.0. 



4. CONCLUSIONS 

We have carefully investigated, based on the latest SNe la data, BAO and SGL statistic, 
the constraint on the equation-of-state parameter w of dark energy for both constant and 
time varying cases in the flat cosmology. The influences of Hubble constant h and matter 
density Qm on the fitting results are carefully demonstrated. The typical redshift measured 
by the three kinds of observations is 2; ~ 1 and far smaller than the redshift of CMB involved, 
their constraints on the parameter w are effective and significant only for the redshift region 
z < 1.5. 

The influence of the equation-of-state parameter w on the density p{z) of dark energy 
in the universe and the distance d{z) makes SNe la data a powerful probe of dark energy. In 
this paper, we have utilized the latest 324 SNe la data provided by Hicken et al. (2009) using 
MLCS17 light curve fitter with the best cuts Ay < 0.5 and A < 0.7 to carefully investigate 
the constraint on the constant w and especially the time- varying w{z) = wq + Waz/{1 + z). 
The influences of Hubble constant h and matter density Qm have also been studied carefully. 
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For the constant w case, we have shown that: (a) the best-fit results for the three correlated 
parameters are found to be (^m? w, h) = (0.358, —1.09, 0.647), it is seen that Qm is somewhat 
large in comparison with Qm = 0.26 of the concordance cosmology provided by WMAP five 
year data(Komatsu et al. 2009); note that using a different parameterization of dark energy, 
Huang et al. (2009) presented a best-fitted Qm = 0.446 from SNe la data, which is even 
larger but still consistent with our result at 95% C.L.; (b) the fitted results of both {flM,w) 
from the SNe la data is very sensitive to the value of Hubble constant h: when h increases, 
the fitted increases and w decreases very rapidly; (c) for parameter Q,m, the SNe la data 
alone can only give a large allowed region 0.0 < < 0.5 at 95% C.L.; (d) for the constant w 
case, the hkelihoods are found to be: w — — O.SSIqJq and — ^-^^-oaI ^^^^ marginalizing 
other parameters in obtaining each of them, which is consistent with the ACDM at 95% 
C.L.; (c) there is clear degeneracy of constant w and Qm, which depresses the power of SNe 
la to constrain two parameters. It has been shown that the parameter w is limited from 
—2.0 to —0.5 and VIm should be less than 0.5. 

In particular, we have paid special attention to the constraints on the time-varying 
case parameterized by two parameters (wo, Wa)- After marginalizing the parameters Qm and 
h, we have obtained the fitting results {wQ,Wa) = (— 0.73^Qgy, 0.84]'^j;Q^34), which indicates 
that (a) the SNe la data alone have only a poor constraint power on the parameter Wa, an 
extra restriction of Qm is necessary, so that the constraint of SNe la on the parameters wq 
and Wa can be much improved; (b)the likelihood of parameter Wa has a high non-Gaussian 
distribution. 

The summary parameter of BAO can provide a standard ruler by which the absolute 
distance of 2; = 0.35 can be determined with 5% accuracy. This ruler can be a prior and 
has been used to constrain dark energy. The strong gravitational lensing (SGL) statistic is 
a useful probe of dark energy. Through comparing the observed number of lenses with the 
theoretical expected result, it enables us to constrain the parameter w. We have used the 
latest SNe la data together with the BAO (and the CLASS statistical sample) to constraint 
dark energy. For the constant w case, the results obtained from (SNe + BAO) and (SNe 
+ BAO + SGL) only have a slight difference: (a) for the (SNe + BAO) case, the best fit 
results of the three parameters (^m, w, h) are (0.287, —0.885, 0.65) and for the (SNe -|- BAO 
-I- SGL) case, the best fit point is {VLm, w) = (0.298, —0.907) where h has been marginahzed; 
(b) the fitting results of the parameter Q.m are found to be Q,m — 0.29lo;o3 the (SNe -|- 
BAO) case and Qm = 0.29lo:o3 for the (SNe + BAO + SGL) case; (c) the fitting results 
for the constant w case are found to be w = — 0.88l'lo;o9 fo^ the (SNe -|- BAO) case and 
w = -OMtoio for the (SNe + BAO + SGL) case, which are consistent with the ACDM at 
95% C.L.; (d) the most allowed region of parameter w is above the line w = —1. Comparing 
with the fitting results from the SNe la data alone, we have found: (a) the allowed region 



- 15 - 



at 95% C.L. for Qm is reduced to one-fifth; (b) the best fit value of w is almost not changed 
but its variance is reduced very much. 

For the time- varying case w{z) after marginalizing {flM,h), we have obtained the fit- 
ting results {wo,Wa) = (-0.951^:1^,0.411^;^^) for the (SNe + BAO) case and {wo,Wa) = 
(-0.921°:^^, 0.351:^:^1) for the (SNe + BAO + SGL) case. It has been seen that the adding 
of the SGL data makes the constraints on parameter (wo,Wa) to be much improved. For 
both cases, the most allowed region of Wa is above Wa = 0, which indicates that the phantom 
type models are disfavored. Comparing with the fitting results from the latest SNe la data 
alone, we have observed that: (a) the best fit values for wq are decreased by over 0.2 and the 
variances are approximately reduced to one-fourth; (b) the best fit values of Wa are decreased 
by 0.49 and the variances are reduced to one-twelfth. 

In conclusion, the latest MLCS17 data set given by Hicken et al. (2009) has provided 
an interesting constraint on the cosmological parameters. The joint analysis of SNe la and 
BAO breaks the degeneracy between w and and leads to a more stringent to constrain 
on the dark energy and matter density than the SNe la data alone, especially for the time- 
varying case with parameters {wo,Wa)- The adding of the SGL data can further improve 
the constraint for {wo,Wa)- A large number of SNe la samples with reduced systematical 
uncertainties in the near future, together with possible new observations on BAO and SGL 
statistic, would be very useful to understand the properties of dark energy and Hubble 
constant. 
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Fig. 1. — The constraint results on the cosmological parameters {^m,w), {h,w) and (/t, from 
top to bottom, respectively. In the upper panel, three pairs of contours are for the given Hubble 
constants h = 0.62, 0.645, and 0.67 from top to bottom, respectively. The crosshairs mark the 
best fit points {nM,w) = (0.01,-0.41), (0.28,-0.87) and (0.48,-2.34) from top to bottom. In 
the middle panel, three pairs of contours are for the given 0,m = 0.35, 0.30, and 0.25 from left to 
right, respectively. The crosshairs mark the best fit points (w,h) = (—1.05,0.647), (—0.91,0.646) 
and (—0.81,0.645) from left to right. In the lower panel, three pairs of contours are for the given 
w = —0.7, —0.9, and —1.1 from left to right, respectively. The crosshairs mark the best fit points 
{nM,h) = (0.19,0.644), (0.29,0.646) and (0.36,0.647) from left to right. 
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Fig. 2. — 68% C.L. and 95% C.L. allowed regions of {^m,w), {h,w) and (/?,, $7m) respectively. 
The crosshairs in three panels markThe best-fit points are {Q.m-,w) = (0.30,-0.92), {w,h) = 
(-0.86,0.645) and (Om,/i) = (0.35,0.647). from top to bottom. 
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Fig. 3. — The likelihoods of the parameters 0,m, w and h. The maximum likelihood points are 
located at JIm = 0.36, w = —0.88 and h = 0.65, respectively. 
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Fig. 4. — The {wo,Wa) contours of SNe la data alone. The crosshairs mark the best-fit point 
{wo,Wa) = (-0.73,0.84). 




Fig. 5. — The {Qm,wo) contours of SNe la data alone. The crosshairs mark the best-fit point 
{nM,wo) = (0.45,-0.68). 
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Fig. 6. — The {CIm, Wa) contours of SNe la data alone. The crosshairs mark the best-fit point 
(J^M,w^a) = (0.44,-4.63). 
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Fig. 7. — The likelihoods of the parameters wq and Wa- The maximum likelihood points are located 
at Wq = —0.8 and Wa = 0.4, respectively. It can be seen that the likelihood of parameter Wa has a 
high non-Gaussian distribution. 




Fig. 8. — 68% C.L. and 95% C.L. allowed regions of {0,m,w) from the joint analysis of (SNe + 
BAO). The best-fit result is {nM,w) = (0.28,-0.88). 
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Fig. 9. — 68% C.L. and 95% C.L. allowed regions of {Qm,w) from the joint analysis of (SNe + 
BAO + SGL) (solid lines) in comparison with the joint analysis of (SNe + BAO) (dotted lines). 
The best-fit result from (SNe+ BAO + SGL) is inM,w) = (0.29, -0.91). 
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Fig. 10. — The likelihoods of the parameters Qm, w and h from the joint analysis of (SNe + BAO) 
and (SNe + BAO + SGL), respectively. The maximum likelihood points are located at JIm = 0.29, 
w = -0.88 and h = 0.65 for (SNe + BAO) and J^m = 0.296 and w = -0.91 for (SNe + BAO + 
SGL). 
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Fig. 11. — the 68% C.L. and 95% C.L. allowed regions of {wo,Wa) from the joint analysis of 
(SNe + BAO) and (SNe + BAO + SGL), respectively. The crosshairs mark the best-fit point 
{wo,Wa) = (-0.95,0.41) for the (SNe + BAO) case and {wo,Wa) = (-0.92,0.35) for the (SNe + 
BAO + SGL) case. 
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Fig. 12. — the likelihoods of parameters wq and Wa from the joint analysis of (SNe + BAO) and 
(SNe + BAO + SGL), respectively. For (SNe + BAO) case, the maximum likelihood points are 
located at wq = —0.94 and Wa = 0.0, respectively. For (SNe + BAO + SGL) case, the maximum 
likelihood points are wq = —0.91 and Wa = 0.34, repectively. 



